gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_2.m

    num=5;
den=conv([1 2],[1 2 1]);
w=logspace(-2,3,100);
[mag,phase,w]=bode(num,den,w);
magdB=20*log10(mag);
subplot(2,1,1);
semilogx(w,magdB);
grid;
title('ϵͳBodeͼ');
xlabel('Frequency(rad/sec)');
ylabel('Magnitude(dB)');
subplot(2,1,2);
semilogx(w,phase);
grid;
xlabel('Frequency(rad/sec)');
ylabel('Phase(deg)');